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Abstract 

We study glueball G production in gluonic penguin decay B XgG, using next-to- leading 
order b sg* gluonic penguin interaction and an effective coupling of a glueball to two gluons. 
The effective coupling allows us to study the decay rate of a glueball to two pseudoscalars in the 
framework of chiral perturbation theory. Identifying the /o(1710) to be a scalar glueball, we then 
determine the effective coupling strength with the branching ratio of /o(1710) — > KK. We find 
that the charm penguin to be important and obtain a sizable branching ratio for Br(S XsG) in 
the range of (0.7 ~ 1.7) x lO"'^. Rare hadronic B decay data from BaBar and Belle can provide 
important information about glueballs. 
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The existence of glueballs is a natural prediction of QCD. However, glueball state has 
not been confirmed experimentally. The prediction for the glueball masses is a difficult task. 
Theoretical calculations indicate that the lowest lying glueball state is a scalar with a mass 
in the range of 1.6 to 2 GeV. Recent quenched lattice calculations give a glueball mass ma 
equals 1710 ± 50 ± 80 MeV [ij. These results support that the scalar meson /o(1710) to 
be a glueball. Phenomenologically, /o(1710) could be an impure glueball since it can be 
contaminated by possible mixings with the quark-antiquark states that have total isospin 
zero li, i, i, H, i, 0, B . These mixing effects can be either small 0, 0, H, or large 0, i, i, B , 
depend largely on the mixing schemes one chose to do the fits and complicate the analysis. 
For simplicity, we will ignore these mixing effects and assume /o(1710) is indeed a glueball 
throughout this work. 

Since the leading Fock space of a glueball G is made up of two gluons, production of 
glueball is therefore most efficient at a gluonic rich environment like J/ip otT (qq )'y —>■ Gj 
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Direct glueball production is also possible at the e~^e~ [11| and hadron [1^ colliders. 
In this work, we point out another interesting mechanism to detect a glueball via the rare 
inclusive process B —>■ XgG decay. The leading contribution for this process is shown in Fig. 
[1], where the squared vertex refers to the gluonic penguin interaction and the round vertex 
stands for an effective coupling between a glueball and the gluons. The gluonic penguin 
b — s> sg* has been studied extensively in the literature and was used in the context for 
inclusive decay b sgrj' [1^. The effective interaction for b — > sg* with next-to- leading 
QCD correction can be written as flT 



Gf gs 
7f 4^ 



V:,Vt^-s{p') [AFi(g2^^ -q^^L- tm.F^a^^q" R] T'^ 



(1) 



where AFi = An^G^^q, fi)+GQ{q, fi))/as{fi) and F2 = —2Gs{f^) with Cj(g, fi) {i = 4, 6, and 8) 
the Wilson's coefficients of the corresponding operators in the AB = 1 effective weak Hamil- 
tonian, q = p — p' = k + k', and T"^ is the generator for the color group. We will use the 
next-to- leading order numerical values of AFi and F2 IJ] . The top quark contribution gives 

charn 
= Pf^'^'^(g,/i) with 



^^top _ _^_gg gj^^ p^"p _ -|_o.288 at yU = 5 GeV; whereas the charm quark contribution 



involves a q dependence through C4 /i) 



C charm 
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(2) 



and 



/' 

^0 



Q{q,m,n) = 4 / dxx{l — x)\n 



— x{l — x)q^ 



2 [3g2 ln(mVu2) - 12m2 - Sg^l Ai2m^ + q^) JArn^ - q^ 

H ^= arctaui 



9g2 



Am? — q^ 
(3) 



Here mc = 1.4 GeV is the charm quark mass and C2(/i = 5 GeV) = 1.150. 

The following effective coupling between a scalar glueball and two gluons was advocated 
recently by Chanowitz [l5| 



£ = fGG^G'^''^ 



(4) 
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where Q is the interpolating field for the glueball G, G^j, is the gluon field strength, and 
/ is an unknown coupling constant. This form of effective coupling suggests that glueball 
couples to the QCD trace anomaly. 




FIG. 1: The leading diagram for B XgG. 



The interaction of a scalar glueball with light hadrons through the trace anomaly can be 
formulated systematically by using techniques of chiral Lagrangian, as was briefly mentioned 
in The kinetic e nerg y and the symmetry breaking mass terms for the light pseudoscalar 
mesons are given by |l7| 

= S.[TT{d^W,J:^) + Tri^X^^ + exO], (5) 



where = 132 MeV being the pion decay constant, 
SU (3) pseudoscalar octet meson. 



S = exp(2zn//^) with U the 
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and 



X = 2-Bo diag(m, m, m^) = diag(m^, m^, 2m^ 



(6) 



(7) 



with Bq = 2031 MeV. Here we have neglected the isospin breaking effects due to small 
mass difference between the light u and d quarks and used rh = rriu = rrid. The QCD trace 
anomaly is well known and given by [13] 



ha. 



in 



(8) 



where 6 = 11 — 2nj/3is the QCD one- loop beta function with nj = 3 being the number of 
light quarks. Treating the effective interaction (jl]) as a perturbation to the energy momentum 
stress tensor, one would then modify to be 



bos 
8^ 



G^G^"" (i + f^g^ + E (^-^ + f'l^^ 1^ 



(9) 
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where fg = /asmq(167r a/2/3/?) In [(1 + — P)] being the one-loop induced Gqq couphng 



15 | with /? = (1 — 4m^/m^)^/^. Note that fg is proportional to /. The corresponding chiral 



Lagrangian is thus modified accordingly 

+ g/l Tr [(Hx + 2So/xe)«* + ax + aSoAS)?] , (10) 

where = diag(/„, fd, fs)- Using the above chiral Lagrangian, one can then calculate the de- 
cay rates for G — > vt+tt^, vr^vr^, K^K~ , K^K , rf^Tf , and 7r°?7°. Since /o(1710) is interpreted as 
the glueball, one can use the experimental branching ratio Br(/o(1710) KK) = 0.38^q^I 
and its total width 137 ± 8 MeV to be our input. Together with the decay rate for- 
mulas derived from the chiral Lagrangian (fTOl) and a value of the strong coupling constant 
as = 0.35 extracted from the experimental data of r decay, one can estimate the unknown 
couphng / = 0.07^°:°;^ GeV-^ 

In our estimation of / given above, we have extrapolated low energy theorems to the 
glueball mass scale. One might overestimate the hadronic matrix elements in due course. 
This implies that the extracted value of / would be too small. 

Decay G —>■ KK has also been obtained using perturbative QCD calculations [l6|. This 
approach can also give some estimate of the amplitude. The problem facing this approach is 
that the energy scale may not be high enough to have the perturbative QCD contribution to 
dominate. Using the asymptotic light-cone wave functions, we find that for a given branching 
ratio for G — > KK decay, the resulting / would be about 20 times larger than the value 
obtained above using the chiral approach. Incidentally, the value of / derived from the chiral 
Lagrangian is within a factor of 2 compared with the value estimated just by using the free 



quark decay rate of G — >■ ss [12|]. In our later calculations, we will use the conservative value 
/ = 0.07!l;og°8 GeV~^ determined using the chiral Lagrangian given in Eq. lfTOj) . This will 
give the most conservative estimate for the branching ratio since the chiral approach gives 
the smallest /. 

With the two effective couplings given in Eqs.([T]) and (jl]), the following decay rate for 
b sgG (Fig.([T])) can be obtained readily 



*sgG 





Glml 




V 4iV, J 
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X dx dy\ I AFil^co + Re(AFiF2*)- + 1^21^4 \ , (11) 

Jo Jv^ L y y ) 



with 



y± = -^ 



'l-x + x')±^{l-x + x'Y - Ax' . (12) 



In Eq.f llll) . Nc is the number of color and Co,i,2 are given by 

Co = ]^[-2x^y +{l-y){y- x'){2x + y- x')] , 
ci = {1 - y){y - x'y , 

C2 = \[2xV -{l-y){y- x){2xy -y + x')] , (13) 



4 



with X = {p' + k)'^/ml, y = {k + fc')^/m^, and x' = rn^/ml. 

Using the value of / determined above from the chiral Lagrangian, we find the branch- 
ing ratio for h sgG to be 4.5 x 10~^(/ ■ GeV/0.07)^ with just the leading top pen- 
guin contribution to AFi is taken into account. The correction from the charm pen- 
guin is nevertheless substantial and should not be neglected. Inclusion of both top and 
charm penguins gives rise to an enhancement about a factor of 3 in the branching ratio 
Br(6 sgG) ^ 1.3 x 10~'^(/ ■ GeV/0.07)^. Since /o(1710) has a large branching ratio into 
KK, the signal of scalar glueball can be identified by looking at the secondary KK invari- 
ant mass. The recoil mass spectrum of can also be used to extract information. The 
distribution of (iBr(6 sgG) /dMx^ as a function of the recoil mass of Xg is plotted in 
Fig.©. 



Recoil Mass Distribution of b ^ sgG 
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FIG. 2: (iBr(i? XgG) / dMxs in unit of 10 ^. In the figure the h quark mass m?, is taken to be 
4.248 GeV, = 0.21 at the b quark scale, tb = 1.674 x lO^^^ and / = 0.07 GeV"!. 

Analogously one can also study inclusive B decays into a pseudoscalar glueball G, whose 
leading effective coupling to the gluon can be parameterized as [19|], 

C = fgG%G''^^ with g;:, = ^e^..^G"^"^ (14) 

and Q denoting the interpolating field for the pseudoscalar glueball. The decay rate for 
this case can be deduced from the scalar glueball one by replacing the coupling / with / 
and the mass mc with the pseudoscalar glueball mass ttlq in Eg. (11 II) . This also reproduces 
previous result obtained in Ref. [l3| for a similar process h — > sgrj' . We therefore expect 
similar distribution and branching ratio for the pseudoscalar glueball production from the 
B decay as in the scalar case that we have studied in this work. 

Recently, BES has observed an enhanced decay in J/ip —>■ jri'mi with a peak around the 
invariant mass of r^Vvr at 1835 MeV [20|. Proposal has been made to interpret this state 
X(1835) to be due to a pseudoscalar glueball j2l|. Taking X(1835) to be a pseudoscalar 
glueball, we would obtain a branching ratio of 3.7 x 10~^(/ ■ GeV/0.07)^ with top penguin 
contribution only, and is enhanced to 1.1 x 10~^(/ ■ GeV/0.07)^ if charm penguin is also 
included. If the coupling / is of the same order of magnitude as /, the branching ratio for 
B XgG is also sizable. Pseudoscalar glueball may also be discovered in rare B decays. 
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Since a pseudoscalar glueball cannot decay into two pseudoscalar mesons, the identification 
of a pseudoscalar glueball necessitates the study of the three-body system TyVvr. This makes 
the analysis more difficult compared with the case of a scalar glueball. 

To conclude, we have studied inclusive production of a scalar glueball in rare B decay 
through the gluonic penguin and an effective glueball-gluon interaction. The branching 
ratio is found to be of the order 10~^. B decays into a pseudoscalar glueball through gluonic 
penguin is also expected to be sizable. Note that we have used a conservative estimate of 
/. Observation of S ^ Xsfo{l710) at a branching ratio of order 10~^ or larger will provide 
an strong indication that /o(1710) is mainly a scalar glueball. With more than 600 millions 
of BB accumulated at Belle and more than 300 millions at BaBar, test of Br (6 — > sgG) at 
the level of 10"'^ is quite feasible. We strongly urge our experimental colleagues to carry out 
such an analysis. 
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